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Abstract 

Consider axisymmetric strong solutions of the incompressible Navier-Stokes equa- 
tions in M.^ with non-trivial swirl. Such solutions are not known to be globally defined, 
but it is shown in [M1[T] that they could only blow up on the axis of symmetry. Let 
z denote the axis of symmetry and r measure the distance to the z-axis. Suppose 
the solution satisfies the pointwise scale invariant bound |t'(a;,f)| < C*(r^ — t)~^/^ for 
—To < t < and < C* < oo allowed to be large, we then prove that v is regular at 
time zero. 

1 Introduction 

The incompressible Navier-Stokes equations in cartesian coordinates are given by 

dtv + {v ■ V)v + Vp = Av, div v = 0. (N-S) 

The velocity field is v{x, t) = (^1,^2,^3) : x [-Tq, 0) and p{x, t) :R^x [-Tq, 0) ^ M 

is the pressure. It is a long standing open question to determine if solutions with large 
smooth initial data of finite energy remain regular for all time. 

In this paper we consider the special class of solutions which are axisymmetric. This 
means, in cylindrical coordinates r,6,z with {xi,X2,X3) = {r cos9,r sinO, z), that the solu- 
tion is of the form 

v{x,t) = Vr{r,z,t)er +vg{r,z,t)ee + Vz{r, z,t)ez. (1.1) 
The components Vr, vg,Vz do not depend upon 9 and the basis vectors er,eo, Cz are 

— ,— ,0 , ee=[ ,— ,0 , e^ = (0,0,l). 

r r / \ r r / 



'National Taiwan University and National Center for Theoretical Sciences, Taipei Office, email: 
chchchen@math.ntu.edu.tw 

^Harvard University, emaif: strainQmath. harvard.edu 
■''University of British Columbia, email: ttsai@math.ubc.ca 
^Harvard University, email: htyau@math.harvard.edu 



1 



The main result of our paper shows that axisymmetric solutions must blow up faster than 
the scale invariant rate which appears in (jl.2p below. 

For R > define B{xq, R) C as the ball of radius R centered at xq. The parabolic 
cylinder is Q{Xo,R) = B{xo,R) x (to - -R^,io) C centered at Xq = {xQ,to). If the 

center is the origin we use the abbreviations Bj^ = B{0, R) and Qji = Q{0, R). 

Theorem 1.1 Let {v,p) be an axisymmetric solution of the Navier-Stokes equations (|N-S|) 
in Z? = M'^ X (— To,0) for which v{x,t) is smooth in x and Holder continuous in t. Suppose 
the pressure satisfies p G L^/^{D) and v is pointwise bounded as 

\v{x,t)\ <C4r^ -t)-^/^, {x,t)eD. (1.2) 

The constant C^, < oo is allowed to be large. Then v G L'^{Br x [— To,0]) for any R> 0. 

We remark that the exponent 5/3 for the norm of p can be replaced. However, it is the 
natural exponent occurring in the existence theory for weak solutions, see e.g. |2gB61 [26]. 

mm- 

Recall the natural scaling of Navier-Stokes equations: If {v,p) is a solution to ()N-Sp . 
then for any A > the following rescaled pair is also a solution: 

v^{x, t) = Xv{Xx, X^t), p^{x, t) = X^p{Xx, XH). (1.3) 

Suppose a solution v{x,t) of the Navier-Stokes equations blows up at Xq = (xo,to). Leray 
|ig[L6l[T6] proved that the blow up rate in time is at least 

\\v{;t)\\L^>C{to-t)-'/\ 

Caffarelli, Kohn, and Nirenberg [MH] showed that for such a blow-up solution the average 
of \v\ over Q{Xq,R) satisfies 

[ \v\^ + \pf/^dxdt) ' > ^. 

.Nr\ Jq{Xo,r) J ^ 

See also [TM9lliyy2y2ll213B7IB7l [37] . Thus, the natural rate for blow-up is at least 



Both this and the rate (|1.2p are invariant under the natural scaling p.3p . 

The Serrin type criteria [2gMl2a[Lai3l[T3gm[gi/i7..7.2S29ll2iJB]]|jlimi5Bl [8] states that v is 
regular if it satisfies 

3 2 

II^IIl|I,|(Qi) < oo, - + - < 1, s,g G (2,oo), or (s, g) = (2, cxd). (1.5) 
q s 

Above, for a domain C M'^, we use the definition 

\H\LILUnx{tut2)) ■= llll^(^'0llL|(f7)llL?(ti,t2)- 
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For any Xq = (xo,to) £ Qii (|l-5p implies the following local smallness of v: 



lim||^llLfL|(Q(Xo,iJ)) =0. (1.6) 

Therefore (jl.Sp is a so-called e-regularity criterion since it implies that the norm is locally 
small. For {q.,s) = (3, oo), (jl.6p does not follow from (jl.5p . Hence the (g, s) = (3, oo) end 
point regularity criterion (jl.Sp proved in |27IB7L[2^M1 [4] is not an e-regularity type theory. 

However these criteria do not rule out blow-up with the natural scaling rate (|1.4p . It 
is a fundamental problem in the study of the incompressible Navier-Stokes equations to 
determine if solutions to (IN-Sp with the following scale invariant bound are regular 

If a self-similar solution satisfies this bound then it is known to be zero 133651 135] (the 
self-similar solution from [21E21 I22] belongs to LfLl). 

Theorem 11.11 rules out singular axisymmetric solutions satisfying the bound (jl.7p . In 
fact ()1.2p is considerably weaker than ()1.7p and is also not a borderline case of the Serrin 
type criterion. For example (|1.2p implies that v £ L'^{Qi) for (7 < 4, but not for q £ [4,5). 
The borderline of the Serrin type criterion, on the other hand, is f S L^{Qi). 

We now recall the previous results on the regularity of axisymmetric solutions to the 
Navier-Stokes equations. Global in time regularity was first proved under the no swirl 
assumption, vq = 0, independently by Ukhovskii-Yudovich [39661136] and Ladyzhenskaya 
|14[L41 [n] . See |igiL51 [T5] for a refined proof and |1]|[LH [TT] for similar results in the half space 
setting. 

When the swirl component vg is not assumed to be trivial, global regularity is unknown. 
But it follows from the partial regularity theory of [IltUfT] that singular points can only lie on 
the axis of symmetry. Any off axis symmetry would imply a whole circle of singular points, 
which contradicts [MKl]. Neustupa-Pokorny |23B3l[23^4B4l[24] proved regularity assuming 
the zero dimensional condition Vr S LfL% with 3/q + 2/s = 1, 3 < q < oo. Regularity 
criteria can also be put on the vorticity field uj = curlv: 

Uj{x,t) = UJrer + i^eeg + LOzCz, (1.8) 

where 

UJ^ = -dzVg, LUe = dzVr - OrVz, UJz = {dr + r~^)vg. 

Chae-Lee |2Bt[2] proved regularity assuming finiteness of another zero-dimensional integral: 
ujg € LfL% with 3/q + 2/s = 2. Jiu-Xin [IQHOl fTO] proved regularity if the sum of the zero- 
dimensional scaled norms jQ^{R~^\ujg\'^ + R~^\vg\'^)dz is sufficiently small for R > small 
enough. Recently, Hou-Li [QBl lQ] constructed a family of global solutions with large initial 
data. 

The main idea of our proof is as follows. The bound (|1.2p ensures that the first blow up 
time is no earlier than t = 0. For t £ (— Tq, 0) we show that the swirl component vg gains a 
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modicum of regularity: For some small a = a(C*) > 0, (jl.2p enables us to conclude that 

\ve{t,r,z)\<Cr''-\ (1.9) 

We prove (jl.9p in Section [3l This estimate breaks the scaling, thereby transforming the 
problem from order one to e-regularity, which is shown to be sufficient in Section [2j 

2 Proof of main theorem 

In this section we prove Theorem 11.11 First we show that our solutions are in fact suitable 
weak solutions. Then we make use of (|1.9|) . to establish our main theorem. 

2.1 Suitable weak solution 

We recah from pM6l[2MSlllllLga9^ [T9] that a suitable weak solution of the Navier-Stokes 
equations in a domain Q C x M is defined to be a pair {v,p) satisfying 

vGLriliQ), VvGL\Q), pGL^^\Q). (2.1) 

Further {v,p) solve ()N-Sp in the sense of distributions and satisfy the local energy inequality: 

2[\Vv\^ip<[{\v\\dtip + Aip) + {\v\'^ + 2p)vV^}, VipGC^iQ), ip>0. (2.2) 

JQ Jq 

To prove interior regularity, we do not need to specify the initial or boundary data. 

We define a solution v{x,t) to be regular at a point Xq if v £ L°°{Q{Xq, R)) for some 
> 0. Otherwise v{x,t) is singular at Xq. We will use the following regularity criterion. 

Lemma 2.1 Suppose that {v,p) is a suitable weak solution of ()N-Sp in Q{Xq,1). Then 
there exists an ei > so that Xq is a regular point if 

limsup^ / \v\^ < ei. (2.3) 

mo -n- Jq{Xo,R) 

This regularity criterion, which is a variant of the criterion in [HKl], was proven in 
[MHlEl]; see [HE1[8] for more general results. The condition (j2.3p does not explicitly involve 
the pressure, but one does require p E L^^'^{Q{Xq, 1)) because the pair {v,p) is assumed to 
be a suitable weak solution. 

2.2 Preliminary estimates 

In this subsection we show that the solution {v,p) in Theorem 11.11 is sufficiently integrable 
to be a suitable weak solution, and we derive estimates depending only upon C* of (jl.2p . 
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We estimate the pressure with weighted singular integral estimates. We therefore first 
estimate v in weighted spaces. Fix /3 G (1,5/3). For t £ {—Tq,0) by (jl.2p we have 

v(x,t)\'^ , ^ f 1 C^rdrdz , , , t ^ ^ 

^ dx< I —— —= I + I + I =h + l2 + h. 



Each of these integrals can be estimated as follows 



|<l,r>l J\z\<l,r<l 



dz f°° C^rdr , , , 

OO ^ 



|/3|< / + ^rdr<c|t|-(^+/^)/2. 
Summing the estimates and using /3 > 1 we get 

[ J_\y(^x,t)\Ux<C+c\t\-^^+^y^. 

Define i?i's to be the Riesz transforms: Ri = J^^ - We consider the singular integral 
p{x,t) = j ^aj(9-,(7;it;j)(?/) ^^^J"_^^ dy = RiRj{viVj). 



To show that this singular integral is well defined for every t, we use the L'?(]R'^) -estimates 
for singular integrals with Ag weight 130601 130] . Specifically, we use q = 2 and the A2 weight 
function We have the estimate 

\j,\p{x,t)\'^ dx < c [ -^\v{x,t)\Ux<c + c\t\-^^+'^^/\ (2.4) 



Choose 7 G (1/2 + 5/3/6, 3). Holder's inequality gives us the bound 



\ 5/6 / \ 1/6 



|a:|>l l^l"^ W|a;|>l \x\^ I \J\x\>l 



< 00. 



We will use these bounds to show that the pressure p can be identified with p. 

Let h{x, t) = p{x, t) —p{x, t). Then h is harmonic in x, ^xh{x, t) = 0, and by assumption 
p{-,t) G L^/'^(R^) for almost every t. For each such t we have 



f |/l(x,t)|5/3 /" , , , f \p(x,t)\^/^ , 

/ ' \Y dx<c \p{x,t)\^^-^dx + c '^^ ' /' dx<QO. 

J\x\>l \xP J\x\>l J\x\>l fP 

We may thus conclude from using a Liouville theorem that h{x,t) = for all x if 7 < 3. 
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To see the last assertion, fix a radial smooth function > supported in 2 < |x| < 4 
satisfying J cp = 1. For any j; G M'^ with R > \x\ we have 

h{x,t) = j h{y,t)R~^4>{x + y/R)dy. 
This is the mean value theorem for harmonic functions. Define A = B^n — Bf{, then 
\hix,t)\ < cR-^ Jjh{y,t)\dy < c/?-3+(6+37)/5 Qjy\-'r\h(y,t)\'/3dyy\ 

This clearly vanishes as i? — > oo. Thus p{x, t) = p(x, t) for all x and for almost every t. 

Next we show that {v,p) form a suitable weak solution. From Holder's inequality, ([27 
and P < 5/3 we conclude that 



/ 



\p{x,t)\^/'^ dxdt < c j (^j dt<c. (2.5) 

The pointwise estimate ()1.2p on v implies 

veLlLliQ^l - + ->^. (2.6) 
q s 1 

We will use (s, q) = (3, 3). We also see from (jl.2p that v E L'^{Bi x (— Tq, — e)) for any small 
e > 0. Thus the vector product of (IN-S|) with rnp for any ip £ C^{Qi) is integrable in Qi 
and we can integrate by parts to get the local energy inequality (|2.2p with Q = Qi. In fact 
we have equality. 

Now, for any R G (0, 1) and £ (— i?^, 0), we can choose a sequence of ip which converges 
a.e. in Qji to H{to — t), the Heviside function that equals 1 for t < to and for t > to- Since 
the limit of dt(p is the negative delta function in t, this gives us the estimate 

vix,t)\^dx+ [ \Vv\^<Cr[ i\vf + \pf/^). (2.7) 



ess sup 

-m<t<QJBR JQr 

These estimates show that {v,p) is a suitable weak solution of (|N-Sp in Qji. Note that these 
bounds depend on C* of (11. 2p only, not on ||p||i5/3(iR3x(_To,o))- 

2.3 Scaling limit 

To show Theorem II. H it suffices to show that every point on the z-axis is regular. Sup- 
pose now a point x^, = (0,0, X3) on the z-axis is a singular point of v. We will derive a 
contradiction. Define X^: = {x^,0). Let (t;'^,^'^) be rescaled solutions of ()N-Sp defined by 

v^{x,t) = XviX{x - x^),X'^t), p^{x,t) = \^p{X{x - x^),\^t). (2.8) 
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By Lemma \2.1\ there is a sequence Afc, k gN, so that Afc ^ as A; ^ oo and 

v^^\^ = i / \v\^ > ei. (2.9) 

We will derive a contradiction to this statement. 

For {v^,p^) with < A < 1, the pointwise estimate (11. 2p is preserved: 

< C,(r2 -t)-i/2^ eM^ x (-ro,0). 

We also have by rescaling 

The argument in the previous subsection provides the uniform bounds for g G (1)4): 

/ b^r + l/l^/^ < esssup f \v^{x,t)\^dx+ f \Vv^\'^<C. (2.10) 
jQi -R?<t<QJBii JQa 

Above the bound for px follows from (|2.5p . the bound for follows from (|1.2p . and the 
energy bound then follows from (|2.7p . 

Thus from the sequence A^ we can extract a subsequence, still denoted by Afc, so that 
{y^k ^p>^k^ weakly converges to some limit function {v,p) 

v^^^v in L^iQn), Vv^'^ ^Vv in L'^{Qr), p^'^ ^ p in L^/^{Qr). 

Moreover since {v^,p^) solves ([N-5]) with bound (l230]) . we also have the uniform bound 

||5jl)^||^3/2((_R2^0);H-2(B^)) < C. 

We can then apply Theorem 2.1 of |33B3( [33l chap. Ill] to conclude that the f^* remain 
in a compact set of L^/^(Qk). Therefore (a further subsequence of) v'^'' — > v strongly in 
L^/'^{Qr). Since the v^^ remain bounded in L'^{Qr) for all g < 4, we deduce that v^^ v 
strongly in L'^{Qr) for all 1 < g < 4. 

2.4 The limit solution 

The convergence established at the end of Section 12.31 is sufficient to conclude that the 
limit function {v,p) is a suitable weak solution of the Navier-Stokes equations in Qr, as in 
p mill 0.911 9l [T9]. Since v satisfies ()1.2p so does v. Hence v is regular at any interior point of 
Qr, and t = is the first time when v{x,t) could develop a singularity. 

To gather more information we use axisymmetry. We will argue in this subsection and 
the next that the estimate (11. 9p (proven in the Section [3]) is enough to conclude that our 
solution is regular. In particular (jl.Op tells us that 

< CA" ^ as A i 0. 



/ 
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Thus the hmit v has no-swirl, vq = 0. 

Let (D = V X be the vorticity of v. The component of (D, u)q = dzVr — dj-Vz, solves 

l\ V 
dt + b-V-A + ^]ue- — CJ0 = 0. 



We have used vq = 0. Above 

b = V = Vj-Cr + VzCz, 6 • V = Vr-dr + Vzdz, div 6 = 0. 
We record the Laplacian for axisymmetric functions 

A-^ -— — 

Next define Vl = uig/r. Then solves 



dt + b-V - A--dr]n = f). (2.11) 



r 



We now derive L"^ estimates on using estimates for the Stokes system. 

Since v satisfies (jl.2p . it also satisfies (j2.6p . We will use both = (5/2,5) and 

{s^q) = (5/4,5/4) below. We rewrite (jN-Sp as a Stokes system with force 

{dt - A)vi + Vip = djfij, div-y = 0, fij = -ViVj. 

By the interior estimates of Stokes system (shown in the Appendix) we have 

Hence Q has the bound 

||f^lL.0/19(Q,/,) < l|V^II^5/4^5/.(Q^^^^||l/r||^^^^.0/n(^^^^^ < C. (2.12) 

In Section [2.61 we obtain O G L°° from ()2.1ip . (I2.12p and a local maximum estimate. Then 
in Section 12.71 we show that this is sufficient to conclude Theorem II. li 



2.5 Energy Estimates 

We derive parabolic De Giorgi type energy estimates for (12. lip . To do this we assume that 

\Kr,z,t)\ <a/r. 

This assumption on b is substantially weaker than the one from Theorem ll.il 
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Consider a test function < ("(x, t) < 1 defined on Qji for which = on dBji x [— i?^, 0] 
and C = 1 on QaR for < a < 1. Define {u)± = maxjitn, 0} for a scalar function u. Multiply 
(IXTTD by p{n - k)^^C for 1 < p < 2 and A; > to obtain 





Br 



4(p - 1) 
P 



\v{{n-krj'0\ 



dc 



(^-^)^{C^ + |vcr + 



-^CAC - 2C— + 6 • CVC 

p r 



dt 



dz c^{n - k)l 



r=0 



Notice that the last term has a good sign. 

To estimate the term involving b we use Young's inequality 



4&C • VC < -5 



1 + e 



vie\b\'+^+cs 



g^-2+(l+e)A 

1 + e 



2 a2 



c 



(i+^)A 



This holds for small 5 > and e > to be chosen. Further choose C to decay like (1 — 
near the boundary of S^. If n is large enough (depending on e) we have 



Cs 



1 + e 



2 a2 



c 



< CR-^ 



Br 



We also use the Holder and Sobolev inequalities to obtain 



1 + e 



C2|6|1+^<5LR(-1+^)3/2 / |^|(l+.)3/2 



\ 2/3 



< / |V(t;±C)l 



2 



|V(i;±C)P 



The last inequality is satisfied for example if |6| < C^:/r and e < 1/3. We conclude 



/ vlbC-VC<5C [ \V{v±C)\^ + CR-^ [ ' 

J^3 Jbr 



(2.13) 



The key point which we used here to control the more singular drift term was to split b from 
the main part of the term v±C, using the Young and Sobolev inequalities instead of standard 
techniques which utilize the Hardy inequality type spectral gap estimate to control 
in one step. We choose 6 sufficiently small in order to absorb this term into the dissipation. 

We have drC/r = dpC,/ p since C, is radial; so that the singularity 1//? is effectively 1/R. 
We thus have 



sup / \{^-k)±\p+ [ \V{n-k) 

-a^R^<t<oJB„RX{t} Jq^r 



< 



c 



(1 - ^)'R' Jq 



P/2|2 
± I 



m-k)±\ 



(2.14) 



Our goal will be to establish to bounds for functions in this energy class. 
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2.6 Local maximum estimate 

The estimates in this section will be proven for a general function u = 0, satisfying (|2.14|) : 
Lemma 2.2 Suppose u = satisfies (j2.14p for 1 < p < 2. Then 

sup U± < C{p, C) f i?-3-2 f _ 
Qr/2 V JQr J 

This estimate can be found in |18[L8l[T8] for p = 2. The proof below is similar and we 
include it so that the proof of Theorem 3.1, which uses Lemma 2.2, is self-contained. Our 
choice of p is made merely because those are the ones we need although others are possible. 

Proof. For > to be determined and a positive integer we define 

feTv = 4 = (1 T 2-^) A', Rn = {1 + 2-^)R/2, pn = 
Rn+1 <Rn = {Rn + Rn+i)/2 < Rn- 

Notice that 

Rn-Rn = {Rn - Rn+i)/2 = (2-^ - 2-^-i)i?/4 = p^. 

Define Qn = Q{Rn) and Qn = Q{Rn) C Qn- Choose a smooth test function Cn satisfying 
C7v = lon(5Ar,C = outside Qn and vanishing on it's spatial boundary, < ("a^ <i 1 and 
IVCa^I < Pn^ in Qn- Further let 

A^iN) ={XeQN-- ±{n - kN+i){X) > 0}. 

And An,± = \A^{N)\. Let v± = Cn{u- kN+ifJ^- 
Holder's inequality gives us 



/ 



|(n - kN+i)±\^ < / \v± 



2 



, , , n/(n+2) 

< / |.;±|2(-+2)/n ^2/(n+2)^ 



Qn 

We will use the following parabolic Sobolev inequality: 

2/n 



/ 

JQ 



|n|2("+2)/" < C(n) I sup [ \un / iVnp. 

R \~R'2<t<oJBRx{t} J Jqr 



See |18I18| [T8 \ Theorem 6.11, p. 112]. We are interested in the form 

/ \ 2/n 

[ |nf/2|2(n+2)/n < f f |v„P/2|2. 

JQr \~B'^<t<oJBRx{t} J Jqr 
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As in the above followed by Young's inequality then followed by (j2.14p we obtain 

N n/(n+2) 
|^^|2(n+2)/n \ 

<C[ sup / / |V^;±P 



< C sup 

\-if2,<t<0^B(iJjv)x{t} 



< C7 sup 

V-i?2^<t<0^B(iJjv)x{t} 



JQn J 

I \{U- kN+l)±r + I |V(U - kM+l)i^A 

JB{RN)x{t} JQm J 



Pn JQn 

<9± I \i^u-kM+l)±y<^ f \{u-kN)±\ 

Pn JQn Pn JQn 



Pn JQn Pn JQn 

Further assume > -^~"~^/q(r) I^±|^- And define 



Y^^K-PR-^-^ [ \{u-kN)±\ 
JQn 



Since are increasing for + or decreasing for — and Qn are decreasing, Yat is decreasing. 
Chebyshev's inequality tells us that 

An,± = \{Qn ■■ ±{u - 4+1) > 0}| = : ±{u - 4) > ±(4+1 - 4)}| 

= \{Qn : ±{n - kj,) > i^/2^+n| < 2^^''+'^ R'^+^Yj^. 

Putting all of this together yields 

n/(n+2) 



JQn+1 \JQn / 

\Pn JQn J ^ 

< (^K^R^+^Yn^ (2P(^+i)i2'^+2y^ 



2/(n+2) 
N,± ' 



2/(n+2) 
2/(n.+2) 



.Pn 



We have thus shown that 

v.- . ^ n( AT\\ ^ 



Yn+1 < CiN)Yj, ' 
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Here C{N) = c, 2^(^+3) 22p(^+i)/("+2). We now choose K as 



'0 



Above the constant Co is chosen to ensure that — > as — > oo. q 
2.7 Regularity of the original solution 

The hmiting solution satisfies (j2.1ip , (|2.12p and ()2.14p . We conclude from Lemma 12.21 
that 

^ e i°"(Q5/i6)- 

We further know that curlu = uieg G (Q^/iq) from the above estimate on since ve = 0. 
Also divv = from the equation. Next v £ -L^L^(Q5/;^g) by (|1.2p . We thus conclude 
Vv E LfL^iQi/^) by Lemma [AH Thus v G L°°(Qi/4) by embedding. 

Now we can deduce regularity of the original solution from the regularity of the limit 
solution. Since v G L°°{Qii^ for R sufficiently small we have 

where ei is the small constant in Lemma |2. II Fix one such i? > 0. Since v'^'' v strongly 
in for k sufficiently large we have 



i?2 

But this is a contradiction to (12. 9p . Thus every point on the z-axis is regular; that is, 
there is a radius Rx^ > so that v € L°°{Q{x^, Rxt))- Since any finite portion of the z-axis 
can be covered by a finite subcover of {Q{x^,,Rx^)}, we have proved Theorem II. 1[ 
The rest of the paper is devoted to proving the key Theorem 13.11 



3 Holder estimate for axisymmetric solutions 

We now move from cartesian to cylindrical coordinates via the standard change of variables 
X = [xi,X2,X'i) = (r cos 0, r sin 0, z). For axisymmetric solutions of the form (jl.ip . the 

Navier-Stokes equations ()N-Sp take the form 



+ 6 . Vv, - ^ + = ( A 
ot r or 



dve , V0Vr ( 1 , 



dv, dp 
— + 6 • Vvz + ^ = 



1 d{rvr) ^ dvz_ _ Q 
r dr dz 
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The vector b is given by 

b = VrCr + VzCz, div 6 = 0. 
The equations of the vorticity u = curlf , decomposed in the form p.8|) . are 

duj [ 1 \ 

+ 6 • — OJrdrUr — UzBzUr = (A n]^r 



at \ 

+ t) • Vwe - 2 — d^ue toe = A ^ t^e, 

Ot J. J. y 

+ b ■ ViOz — ^zdzUz — LOrdrUz = Aw^. 



dt 



Although we do not use them. We are interested in the equation for vg, which is independent 
of the pressure. 

Consider the change of variable T = rvg, which is well known (see the references in the 
introduction). The function T is smooth and satisfies 

^ + 6.vr-Ar + -^ = o. (3.1) 

ot r or 

Note that the sign of the term ^ ^ is opposite to that of (j2.1ip . It follows directly from ()1.2p 
that ||r||2,c» < C^; see PE|[2] for related estimates. Since v is smooth, we have T{t,Q,z) = 
for t < 0. The smoothness and axisymmetry assumptions also imply vg{t,0,z) = 0, but we 
will not use this fact. The main result of this section is the following. 

Theorem 3.1 Suppose that T{x, t) is a smooth bounded solution of (13. ip in Q2 with smooth 
b{x,t), both may depend on 9, and 

r|r=o = 0, div6 = 0, \b\<C^/r in Q2. 

Then there exist constants C and a > which depend only upon C^, such that 

\T{x,t)\ < Cllrllioo (Q2)r" in Qi. 



We remark that the condition above is substantially weaker than (jl.2p . and we do not 
need F to be axisymmetric. In the rest of this section, we will prove the theorem. Here 
we are facing two difficulties: First, the condition F|^^q = precludes a direct lower 
bound on the fundamental solution and a Harnack inequality on F (since, when 6 = 0, 
F = is a nonnegative solution which does not satisfy the usual Harnack inequality.) 
Second, the condition 6 < C/r is weaker than the standard assumption 6 < C/|a;| (see 
the discussion below). It turns out that one can develop new techniques incorporating the 
methods introduced by De Giorgi [3E1[3] and Moser |2QB0l[20] to over come these two points. 
However, we do not know if one can follow the approach of Nash [2IIB11I2IIHE1 I6] which relies 
critically on a Gaussian lower bound of the fundamental solution. The proof of Theorem 
3.1 is independent of the rest of the paper. 
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The following related equation has been previously studied by Zhang 139691 139]: 



— + b-Vu-Au = 0. 
ot 

He has shown among other things Holder continuity of solutions to this equation if 6 = b{x) 
is independent of time and b satisfies an integral condition which is fulfilled if say b is 
controlled by His proof makes use of Moser iteration and Gaussian bounds. 



3.1 Notation, Reformulation, and Energy inequalities 

Let X = {x,t). Define the modified parabolic cylinder at the origin 

Q(i?, t) = {X : \x\ < R, -tR^ <t<Q}. 
Here > and r G (0, 1]. We sometimes for brevity write Qr = Q{R) = Q{R, !)• Let 

m2 = inf r, M2 = sup L, M = M2 - m2 > 0. 

Qm Q{2R) 

Notice that 1112 < < M2 since r|j.=o = 0. 

Now we reformulate the problem in Q(2R) into a new function, u, which will be zero 
when |r| is at its maximum value. Specifically, we define 

(2{T-m2)/M if -m2>M2, 
u=s, , , (3.2) 

[2(M2-r)/M else. ^ ^ 

In either case u solves (|3.1|) and < n < 2 in Q{2R). We will further use 

2 / \ 2 

a = u\r=o = — sup |r| = — max{M2, -m2} > 1, 
^ \Q{2R) J ^ 

which follows from our conditions. 

We now derive energy estimates for ()3.ip . Define v± = {u — k)± with A; > 0. We have 
< {2 — k)^ and V- < k. Consider a radial test function < Ci^^t) ^ 1 for which C = 

on dBji X [—tR"^, 0] and ^ < 0. We multiply (j3.ip for u — k with C'^v± and integrate over 

M'^ X [to, i] to obtain 



J to J to 

+ 2TT[{a-k)±f [ [ dz C\=o. 
Jto JR 



14 



We need to estimate all the terms in parenthesis. 

Choose a G (1/4, 1), we require that the test function satisfies C = 1 on Q{aR, r). If we 
further choose C{x,to) = then, using (|2.13|) . we estimate (|3.3|) as follows 



)Jb(ct 



sup / 4 + / |Vf±p 



-To-2R2<t<oyB(crR)x{t} JQ{aR,T) 

< ,n ^**)2m I vl + CrR'iia-k)^] 
r(l - a)^K- JQiR,r) 

If we alternatively choose ^ = C{x) then (jS.Sp takes the form 

sup /" v'i + [ [ \Vv±\^ - [ u| 

to<s<tJB{aR)x{s} JtoJB{aR) J Buy: {to} 



< c 



f f vi+CTR\a-k)±\' 

J to J Bfi 



(3.4) 



(3.5) 



(l-<T)2i?2 

Notice that there is no appearing in this energy inequality (j3.5|) compared to 

The energy estimates (|3.4p and (|3.5p are the standard parabolic De Giorgi classes except 
for the last term. Our goal will be to use them to show that the set where T is very close 
to its largest absolute value or, equivalently, the set where u is almost zero is as small as 
you wish. We establish this fact in the following series of Lemma's. 



3.2 Initial Estimates 

Later on we will use the two standard Lemma's below in a non-standard iteration scheme 
of sorts to show that the set where u is almost zero has very small Lebesgue measure. 

Lemma 3.2 Suppose there exists a to £ K > and 7 G (0, 1) so that 

\{xeBR:u{xM)<K}\ <-i\Br\. 

Further suppose that u satisfies (j3.5p for V- . Then for all 7] £ (0,1 — y/j) and ^ £ (7/(1 — 
r/)2, 1) there exists 6 G (0, 1) such that 

\{x(^BR:u{x,t)<r]K}\ < ^i\Br\, Vt G [t(,M + {t ^e)R\ 

Here 9 depends only on the constants in (13. 5p and 7. 

We note that the proof shows that ^(7) — > as 7 | 1, but if r is sufficiently small then 
we may take = t when 7 is close enough to zero. And if 7 is small, then can be taken 
almost as small. 
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Proof. We consider v- = (u — K)^ and assume without loss of generality that K <1. 
The energy inequality ()3.5p for this function is 



JB{aR)x{t} JBRx{to} [L-a)K JtoJBR 



We have used (a — K)^ = 0. The Chebyshev inequality tells us that 
A=\{xe B{aR) : u{x, t) < i]K} \ ■ {K - r]Kf < [ 

JB{aR)x{t} 

The quantity 1 — is an upper bound for the measure of A'^, which grants the following 
general inequality 

I {x G Br, u{x, t) < T]K} I I {x G B{(jR),u{x, t) < t]K} \ 3 

- r5~\ h (1 - cr ), 



\Br\ - \Br\ 



Now let a be so close to one that (^i2^y2 + (1 — cr'^) < A*- Then, with r fixed, choose 9 small 
enough that the whole thing is < ;U. □ 
The Lemma above shows continuity in time of the Lebesgue measure of the set where u 
is small and the lemma below shows that if the set where u is small is less than the whole 
set, then the set where u is even smaller can be made tiny. This is an extremely weak way 
to measure diffusion. 

Lemma 3.3 Suppose that u{x,t) satisfies (13. 4p forv-.. In addition 

\{x£BR:u{x,t)<K}\<j\BR\, Vt G [to,to + OR^] = I, 

where K,6 > 0, ^ £ (0,1) and Br x I C Q{R,t). Then for all e G (0,1) there exists a 
6 G (0, 1) such that 

\{X gBrxI : u{X) < 6}\ < e\BR x I\. 
Proof. We denote, for n = 0, 1, 2, 3, . . ., 

Anit) = {x£ Br : u{x, t) < 2-"i^} , An = UtelAnit). 

Clearly \An+i\ < |^„| < \Ao\ < -f\BR x /|. And 

\A'n{t)\ = \{xeBR: u{x,t) > 2-"A'}| = \Br\ - |A„(t)| > (1 - j)\Br\. 
Since 7 < 1, we know that A'f^{t) does not have measure zero. 
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We invoke the following well known version of the Poincare inequality. For any v G 

W^'^{Bii) and for any a, /3 E M with a < /3 we have 

|{x £ Br : v{x) < all < , r ^^^^^^^ ~ / |V7;|, 

where C > only depends on the dimension. Let /? = 2~^K and a = 2~"~^A'. We have 

< T / |V«| = . / |V(n - a)-|. 

We use the Cauchy-Schwartz inequality to bound this integral as 



r C2'^~^ R f 

K«l=/|A„„wi<^^y^__^_jv(„-„,-| 

-\A„ - A„+i\''^ ( f |V(u-o 



C2"+iii 

< 



K(l-7)' 

The energy inequality (|3.4|) . with cri? and R replaced by R and 2R results in 
K{l-l) \rR^ Jq(2R,t) J 

Square both sides of this inequality and dividing by \Br x to obtain 

lA+iP ^ C / \An\ _ |A„+i| 

\Br X /|2 - ei^2(l _ \^ \Br X /| |5r X I\ 

Summing in n, we get 

\BrxI\^ -fT{\BRxI? - eK\l-^)^^^\\BRxI\ \BrxI\ 

C ( \Ao\ \An 



eK^{l--if\\BRxI\ \BrxI\ 
< g l-4o| ^ C7 



We complete the proof by choosing n sufficiently large. □ 
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3.3 Estimate on the measure of the set where u is small 

The next lemma allows us to apply all the machinery above. 

Lemma 3.4 There exists a k £ (0, 1) such that < A < min{KT, 1/8} implies 
\{X £ Q{R,t) : u{X) < X^}\ < (1 -4A)|Q(i?,r)|. 
Proof. We establish a contradiction using energy estimates. Suppose the opposite 
\{X G QiR,T) : u{X) < A2}| > {1-4X)\Q{R,t)\. 

Or equivalently 

\{X G Q{R,t) : u{X) > A^ll < 4A|Q(i?,r)|. (3.6) 

This condition will imply a contradiction to the size condition on a > 1. 

We will test the equation (|3.ip with pnP"^C^ for < p < 1 and C ^ 0- Since u = 
sometimes, in general we should test ()3.ip for n + e with p{u + e)^~^C^ and then send e J, 
to obtain our estimates. However, since the result is the same, to simplify the presentation 
we will omit these details. We have 



„_1 .2^^^ 

QiR,r) 







7 cv 









ot 



[ pvP-\\-/\u) =li^Jl / |VK/2^)|2 
Jq{r,t) p Jq{b.,t) 

+ [ 2uP [-|VCP + ^CAC 
JQ{B,t) I P 

[ puP-\H-Vu = - [ 2uPb-CVC = l5, 
Jq(r,t) Jq(r,t) 



'Q{R,t) JQ{R,r) 
[ puP"^C^-drU = - [ Au^CCp/p- [ dt [ dz2(CV)|^=o 

Jq{r,t) f Jq{r,t) J-tr? Jr 

= h + It- 

In the computation of /g we have used Cr/^ = Cp/p^ which follows if ^ = CiP^t) where 
p = \x\ = \/r^~+~z^. Notice that Ylj=i — 0- arbitrary p G (0, 1), we see that I3 and 
It are both non-positive. 

We choose ( = Ci{p)C2it) where Ciip) = 1 iii -B(i?/2) and Ci(p) has compact support 
in Bji; also C2{t) = 1 if t G [— |ri?^, — |r-R^] and C2{t) has compact support in {—tR'^,0). 
Thus Ii = and we have 

We estimate each of the terms I2 through Iq to obtain a contradiction. 
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By the argument in ()2.13p . we have 

< I |V(u-/\)P + ^ / W. 

P JQ(R,r) ^ JQ{R,t) 

Also note VC = in B{R/2) and so the singularity 1/p is effectively 1/R. Thus, 

jQiR,r) fr", ' JQ{R,r) 

Assuming (13. 6p and using < n < 2, we have 

< - 7^ {A2^|Q(i?,r)| +2P(4A|Q(ii,r)|)} < ^{X'^ + X). 

Here C2 = C2(C*). Take p = 1/2 and k = to get < 1, a contradiction. q 

Lemma 13.41 is the starting point of our iteration scheme. From this Lemma we know 
that there is a G [-tR"^, —2XtR'^] so that 

\{xeBR: u{x,h) < A^} I < (1 - 2X)\Br\. (3.7) 

Then apply Lemma 13.21 with i^T = to p.Th to see, for say r] = X and fi = 1 — X, that 

I {x G Br : u{x,t) < A^} I < {1-X)\Br\, Vt G [ti, ti + ^^ii^] = 



Here 0* = 9 A t and ^ is the constant chosen in Lemma 13.21 From here Lemma 13.31 allows 
us to conclude 

\{X€BrxI,: u{X) <6,}\< ^\Br x /,|, 

where e* > is as small as you want and 5* = (5=i,(e*). 

Then, as in (j3.7p . there exists a t2 € (so that t2 < —XtR?) such that 

\{x(^Br: u{xM) < 5*} I < e,\BR\. (3.8) 

Up till now all the small parameters that we have chosen depend upon r. But above e,,, can 
be taken arbitrarily small independent of the size of r. This is the key point that enables 
us to proceed. It is the reason why we are required to do this procedure twice. 

Now suppose 1 — (T^ = 1/4 and choose first r < 1/8 so that C=n,r/(1 — cr)^ < 1/4. Then 
take from ()3.8|) with < 1/16 above playing the role of 7 in Lemma 13.21 Also rj < 1/2. 
With all this, from Lemma 13.21 we can choose /x < 1 so that 

\{x G Br : u{x,t) < r]6^}\ < ti\BR\, Vi G 1*2,^2 + Ti?2] = /. 

Further, it is safe to assume that 6^ < A; we see that t2 < —XtR^ and so [—XtR'^,0] C /. 
Finally apply Lemma 13.31 again to obtain 

\{X G QiR, At) : u(X) < 5}\ < e\Q{R, Ar)|, (3.9) 
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with e > arbitrarily small. This is a key step in what follows. 

Let U = 6 — u, where 5 is the constant from (j3.9p . U is clearly a solution of (j3.1|) and 
U\r=o = 6-a<0. We apply dMD to U on (5(2(i) (with r = 1) to get 



sup / \{U-k)+\^+ [ |V(C/-A;)+|2 



(1 - ayd^ Jg^d) 

This holds for all /c > and cr G (0, 1). So we can apply Lemma 12.21 to conclude 

/ \ 1/2 

sup(5-n)< \iS-u)+\A . (3.10) 

This inequality combined with ()3.9p will produce a lower bound. 

3.4 Regularity from a lower bound 

Let d = y/XrR so that Q{d) C Q{R, Ar). By (fXTU]) and (I33D, 

/ \ 1/2 

V \Qid)\ J 



which is less than | if e is chosen sufficiently small. We conclude 

inf u>—. 

Qid/2) 2 

This is the lower bound we seek. From it we will deduce an oscillation estimate. 
This entails a bit of algebra. We define 

rud = inf r, Md = sup F. 

QW2) Q(d/2) 

Then from (13.21) we have 



2{md - m2)/M if - ms > Ms, 
mf li = < , 
Q{d/2) y2{M2-Md)/M else, 

Notice that both expressions above are non-negative in any case; thus we can add them 
together to observe that 

^<^{M-osc(F,d/2)} 
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Here osc(r, d/2) = — and osc(r, 2R) = M2 — m2 = M. We rearrange the above 

osc(r,(i/2) < (^l-^ osc{T,2R). 
This is enough to produce the desired Holder continuity via the following. 
3.5 Iteration Argument 

Suppose we have a non-decreasing function uj on an interval (0, Rq] which satisfies 

oj{tR) < juj{R), 
with < 7, r < 1. Then for R < Rq we have 

c^(i?)<^(|^)''^(i?o), (3.11) 

where a = log 7/ log r > 0. 

Iterating, as in ()3.11|) . we get, for Cr = (l — |) supQ^;^) F, that 

osc(r, R) < CrR", Vfl G (0, 1), (3.12) 

for a = 21og(l — |)/ log(Ar/16) > 0. Thus T is Holder continuous near the origin. We have 
proved Theorem 13. 1[ 

Appendix 

Here we collect some estimates needed for Section 2. 

Lemma A.l Let B^^ C Br^ C M'^ be concentric with < R2 < Ri- Let v be a vector field 
defined in Br^ . Let 1 < g < 00 and < a < 1. Then for A; = 0, 1, . . . there is a constant c 
depending on R2, Ri, q, a, k so that 

\N''^^u\\li(Bu^) < c||V''divn||i5(B^j + c|| V'' curlu||j^g(B^j + c||u||ii(B^j. 

and 

||V''+^n||c<.(B^^) < c||V''divu||c-(B^j +c||V''curln||c-(Bflj + c||u||li(b^j. 
This is well-known, see [2^122] . 

Lemma A. 2 (Interior estimates for Stokes system) Fix R G (0, 1). Let 1 < s,q < 00 

and f = (fij) G LfL%(Qi). Assume that v G L*L^((5i) is a weak solution of the Stokes 
system 

dtVi - Avi + dip = djfij, divt; = in Qi. 
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Then v satisfies, for some constant c = c{q, s,R), 

W'^vhlLUQn) < cWfhlLliQ,) + ^vhfLKQ^)- (A.l) 
// instead v is a weak solution of 

dtVi - Avi + dip = Qi, div f = in Qi, 

then 

W^MLtLUQR) ^ 49\\LtL%iQ,) + cIHlslHq^)- (A.2) 

An important feature of these estimates is that a bound of the pressure p is not needed in 
the right side. A similar estimate for the time-independent Stokes system appeared in [32B21 
132] . Note that these estimates improve the spatial regularity only. One cannot improve the 
temporal regularity, in view of Serrin's example of a solution v{x,t) = /(t)V/i(x) where 
h{x) is harmonic. 

Proof. Denote by P the Helmholtz projection in M^, {Pg)i = Qi — RiRkSk, where 
Ri is the i-th Riesz transform. Let r = i?^/^ e {R, 1) and choose C{x,t) E C°°(M^), ( > 0, 
( = 1 on Qt and = on x (— oo, 0] — Qi- For a fixed i, define 

Vi{x,t) = j T{x -y,t- s) dj{Fij){y,s) dyds, 

where F is the heat kernel and Fij = fijC, — RiRkifkjC)- The function Vi satisfies 

{dt - A)di = djFij = [Pd,C{fkj)l=i]i, div v = 0. 

The L^-estimates for the parabolic version of singular integrals and potentials (see [171 P. 7^ 
[IIEai5l[2MEni[3Q], also see [33581 138]. [m2l [T2] and their references), and the usual version 
of L'^-estimates for singular integrals, give 

\Nv\\lil%{q^) + MlilKq,) < 4F\\LtLl < cWfhtLUQi)- (^-3) 
Furthermore, for some function p{x,t), 

idt-A)v + Vp = dj{cfij), dWv = o. 

The differences u = v — v and tt = p — p satisfy the homogeneous Stokes system 

dtu — An + Vvr = 0, diw = in Qr- 

Its vorticity uj = curlu satisfies the heat equation {dt — A)uj = 0. Let W = Ct^^ where 
Qr{x,t) = (^(x/r, t/r^). It satisfies 

{dt - A)W = G:= W{dt - A)Cr - 2{dmQr)draUJ. 
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And thus, for {x,t) S Qr^, 



rt r rt 

iOi{x,t) = Wi{x,t) = I j T{x - y,t - s)Gi{y,s)dyds = / / N'My, s) Uj{y, s) dyds 



where, using tOi = —dij^dkUj, 

Hl;i(.y,s) = dyj,,k{nx -y,t- s){dt - A)Cr + 2div[r(x -y,t- s)VCr]} . 
The functions Hx^y are smooth with uniform L°°-bound for {x,t) £ Q^i. Thus 

||curl'u||i^(Q^3) < C\\u\\LHQ^y 

Since divtt = 0, we have for any q < oo, using Lemma |A.1|, 

W'^^W MLUQr) ^ cMLfLUQ^) < c\\v\\lilI(Q) + c\\v\\LtLUQ)- (A.4) 

The sum of ()A.3P and ()A.4p gives ()A.ip . The proof of (|A.2p is similar: one defines 

Vi{x,t) = J T{x - y,t - s) Fi{y,s)dyds, Fi = giC - RiRkigkC) 
and obtains \\"^'^v\\LtLUQi) + II^IIl?lS{Qi) ^ < c\\9\\l''Li(Qi)- One then estimates 

- v)\\liluqr) ™ s^™^ ^^y- □ 
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